Electrically tunable transverse magnetic focusing in graphene 
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Electrons in a periodic lattice can propagate 
without scattering for macroscopic distances de- 
spite the presence of the non-uniform Coulomb 
potential due to the nuclei [Ij. Such ballistic mo- 
tion of electrons allows the use of a transverse 
magnetic field to focus electrons [2j. This phe- 
nomenon, known as transverse magnetic focusing 
(TMF), has been used to study the Fermi sur- 
face of metals [3j and semiconductor heterostruc- 
tures [4j, as well as to investigate Andreev re- 
flection [3j, spin-orbit interaction [5j, and to de- 
tect composite fermions [6l[7]' Here we report on 
the experimental observation of transverse mag- 
netic focusing in high mobility mono-, bi-, and 
tri-layer graphene devices. The ability to tune 
the graphene carrier density enables us for the 
first time to investigate TMF continuously from 
the hole to the electron regime and analyze the 
resulting "focusing fan". Moreover, by applying 
a transverse electric field to tri-layer graphene, 
we use TMF as a ballistic electron spectroscopy 
method to investigate controlled changes in the 
electronic structure of a material. Finally, we 
demonstrate that TMF survives in graphene up 
to 300 K, by far the highest temperature re- 
ported for any system, opening the door to novel 
room temperature applications based on electron- 
optics. 

The concept of TMF can be illustrated by consider- 
ing electrons entering a two-dimensional system through 
a narrow injector (origin in Fig. la). In the presence 
of a magnetic field 5, electrons will undergo cyclotron 
motion with radius Tc and get focused on the caustic (a 
quarter of a circle with radius 2rc) on which the electron 
density becomes singular (Fig. [l^, top). Moreover, the 
specular refiection of electrons at the boundary of the 
two-dimensional system causes a skipping orbit motion 
which results in focal points at integer multiples of 2rc 
along the x-axis (Fig. [l^, bottom). This basic behav- 
ior still holds for electron motion in a solid as long as 
the Fermi surface has cylindrical symmetry [3]. Hence, 
the magnetic field, 5/, required to focus electrons at a 
distance L is 

r^j) f 2hkY\ / 2h\/7rri\ , , 

where p — 1 is the number of refiections off the edge of 
the system (e.g. p = 1 corresponds to direct injector to 



collector trajectory, without refiections), h is the reduced 
Planck's constant, e is the elementary charge, /cp is the 
Fermi momentum, and where we have used /cp = V^Tm, 
n being the carrier density. 

In order to study TMF in graphene, we fabricate 
Hall bar devices based on high mobility mono- (MLG), 
bi- (BLG), and tri-layer (TLG) graphene on hexagonal 
boron nitride (hBN) substrates [8 (see Methods and 
Fig. [l]^). The multi-terminal geometry required to study 
TMF imposes a minimum mean free path of the order of 
several hundred nm, which has only been possible with 
the advent of G on hBN devices [SHTO]. Figure [l]i shows 
the resistivity of a MLG device as a function of density at 
zero magnetic field. The device exhibits a narrow Dirac 
peak with a strong temperature dependence, which in- 
dicates low disorder [ll|il2j. Its field effect mobility is 
~100,000 cm^V~^s~^ at low temperature, correspond- 
ing to a mean free path of ~1 /im. A similar behavior 
is observed for BLG and TLG devices. The high mobil- 
ity and low disorder enable us to probe TMF in these 
devices. 

We employ the measurement configuration shown in 
Fig. [T]^ to probe the focusing of electrons. Current I\ 
is injected through contact i while contact g is grounded 
and voltage is measured at the collector (contact c) rel- 
ative to contact f. The magnetic field is applied normal 
to graphene. Figure [l^ shows the normalized V^iB^n) 
in MLG, at 5 K. Two sets of features are immediately 
apparent: for \B\ > 2.5 T, we observe Shubnikov-de 
Hass oscillations (SdHOs), forming a usual Landau fan, 
as expected from the measurement setup, which is topo- 
logically equivalent to a longitudinal resistance measure- 
ment. While the SdHOs are very pronounced in quad- 
rants 2 and 4 (top-left and bottom-right), they are nearly 
invisible in quadrants 1 and 3, due to the interference 
of different trajectories of electrons propagating coher- 
ently to the collector [H [T3HT5] (see supplementary infor- 
mation). 

In the low field regime, \B\ < 2.5 T, we observe three 
unusual peaks which do not resemble SdHOs. For posi- 
tive density, these peaks appear on the positive B side. 
The location of these features in the B — n plane (see 
Eq. 1) indicates that these peaks can be associated with 
TMF. The peaks arise when electrons are focused onto 
the collector, resulting in a build up of Vc- The first 
peak corresponds to electrons propagating directly from 
the injector to the collector while for the higher order 



FIG. 1: Transverse magnetic focusing schematics, a, Classical trajectories of electrons injected isotropically from the 
origin at B = Bf (top) and B' = 2Bi (bottom, including one bounce off the edge). Electrons are focused at an integer multiple 
of 2rc along the x-axis. b, Cartoons depicting the band structure of MLG at positive (left) and negative (center) electron 
density. Electron's trajectories at a finite B are shown on the right for Fermi energy > (orange line) and Ef < (blue 
line), c, False color atomic force microscopy (AFM) image of a TMF device. In the TMF measurement, contact i injects 
current li into graphene and the voltage Vc is measured at the collector (contact c) relative to contact f. L is the measured 
distance between the centers of contacts i and c. d. Resistance versus gate voltage of MLG at various temperatures measured 
in the usual 4-probe Hall bar geometry, e, TMF spectrum in MLG at 5 K. TMF peaks from first, second, and third modes can 
be observed clearly for \B\ < 2.5 T. The top insets show representative trajectories for each corresponding mode. At higher B, 
SdHOs are also present. 



peaks electrons reflect off the edge before reaching the 
collector (Fig. [l^, top insets). For negative electrons 
propagate away from the collector and hence no focusing 
peak is observed (Fig. [l^, top left inset). As we tune to 
negative density {E^ < 0), the sign of the charge carri- 
ers flips, and therefore B has to be reversed in order for 
the carriers to be focused at the collector. The ability 
to tune density in graphene enables us to investigate the 
dependence of the focusing fields, or "focusing fan" , 
continuously from the electron to hole regimes in a single 
device over a broad density range, which was never done 
in other systems. 

The values of Bf can be readily calculated, since both 
n and L can be obtained from Hall measurements and the 
AFM image of the device, respectively. Figure [2^ shows a 
zoom-in plot of Fig. [l^, where we have superimposed the 
calculated focusing fields (dashed lines) using the mea- 
sured Lmlg = 500 nm. A discrepancy between the cal- 
culated values and the measured peak locations is clearly 
present. Moreover, we find that the observed b\^^ /p de- 
creases as p increases (Fig. [2]i). The finite width of our 
injector and collector (~100 nm in MLG and BLG and 
~240 nm in TLG, see below) could introduce an error in 



the determination of L and subsequently B^. However, 
a more plausible explanation is the effect of charge ac- 
cumulation near the edges owing to the finite size of our 
graphene devices [16 . We find that charge accumulation 
reduces B^ by the same order of magnitude as that re- 
quired to correct for the discrepancy and, in addition, it 
also explains the decreasing b\^^ /p because, for higher 
the carrier's trajectory is closer to the edge, which fur- 
ther reduces b\^^ (see supplementary information). We 
also note that density fluctuations and small-angle scat- 
tering [14 due to impurities could also affect the car- 
rier's path and its focus. However, a lack of knowledge 
of the detailed disorder potential landscape prevents us 
from determining the change in the value of the focusing 
fields. 

We have observed multiple focusing peaks in all of our 
devices, including BLG and TLG (see below), which indi- 
cates that a significant fraction of the electrons get spec- 
ularly reflected off the graphene edge. From the peak 
amplitudes, we can calculate the measured specularity, 
the ratio between the amplitude of the second mode to 
that of the first mode, which offers information on the 
specular reflection of electrons at the graphene edges. 



3 



VJ I. (Q) 
a -200 200 400 



nr-Nnr 



3 

_ 2- 

I 1- 
o 



I -1 

CD 
Q 

-2 
-3 



▼ 


ij 







-1 1 

Magnetic field (T) 



n = -2.79 X 10^2 cm 



100 



-1 1 

Magnetic field (T) 





3- 




2- 


E 




o 


1 - 


o 


>. 




"w 


0- 


c 

CD 




Q 






-1 - 




-2- 



-100 



100 





Ik 


/ f J 





-1 



1 




Magnetic field (T) 
n= -2.78x1 0^2 cm 2 q R 2R 3R 









I: 




2 -1 


c 


) 


1 





-200 -100 100 




1 
Magnetic field (T) 



n=7.51 X 10^2 cm 2 



R 2a 3R 




Magnetic field (T) 



-1 1 

Magnetic field (T) 



FIG. 2: transverse magnetic focusing in MLG, BLG, and TLG at 5 K. a-c, The TMF spectra as a function of density 
and magnetic field for MLG, BLG, and TLG respectively. The dashed lines are calculated focusing fields using Lmlg = 500, 
i^BLG = 775, and Ltlg = 950 nm, determined from AFM images, d-f, Onsager reciprocal relation in MLG, BLG, and TLG 
respectively. The red and black traces are measured with current and voltage contacts switched (Fig. |2]l, insets). In these 
figures, Bi is the observed focusing field of the first mode. 



We find that the value of specularity ranges from 0.2 
to 0.5 (see supplementary information). It is worth not- 
ing that specularity measurements in semiconductor het- 
erostructures have shown values less than 1 for focused- 
ion-beam etched devices which is similar to our oxygen- 
plasma-etched graphene devices but greater than 1 for 
electrostatically-defined edges [31 [4]. 

We now turn to TMF in BLG. Figures |2]3 and e dis- 
play TMF fans for BLG at 5 K. Evidently, the TMF 
spectra of MLG and BLG are very similar, even though 
their band structures are different (Fig. [2|i-b, insets). 
The similarity arises from the fact that, when only the 
nearest intra-layer 70 and inter-layer 71 hopping param- 
eters are considered, both MLG and BLG have circular 
Fermi surfaces, resulting in the same circular orbit and 
y^-dependence of /cp. The dashed lines in Fig. are 
focusing fields calculated from Eq. 1. For this device, the 
calculated Blf^ for the first order peak are in good agree- 
ment with the measured data, but higher order modes 
show a discrepancy, similar to the MLG case. An addi- 
tional possible source of mismatch in BLG, which does 
not exist in MLG, is trigonal warping [17 of the Fermi 
surface due to the next nearest neighbor inter-layer hop- 
ping term 73. This term transforms the BLG circular 
Fermi surface into a partly triangular surface, altering 



therefore the carrier's trajectory. Hence, in principle the 
values of Bi now depend on the crystallographic orienta- 
tion with respect to the sample axis, and can vary by a 
few tens of mT (see supplementary information). 

Even though TMF cannot be used to differentiate 
MLG from BLG, the TMF spectrum of TLG is remark- 
ably different because of the multiband character of its 
band structure. Figure [2]3 and f show TMF spectra of 
TLG at 5 K, measured at zero electric displacement field. 
Taking only 70 and 71 into account, the band structure of 
TLG consists of a massless MLG-like and a massive BLG- 
like subband at low energy (Fig.[2]3, inset) pl ITMST] . In a 
magnetic field, both subbands give rise to their own TMF 
spectra, with the BLG-like subband having a larger B^^^ 
due to its larger Fermi momentum (for a given E^) • This 
allows us to identify the subband corresponding to each 
peak observed in the data. At high density, the peak from 
the MLG-like subband can be seen at -250 mT (Fig.|2]:3, 
orange dashed line, and small sharp peaks at low field 
in Fig. |2]p) while the peaks from the BLG-like subband 
are visible from —250 mT onward (Fig. [2]3, blue dashed 
lines). We do not observe higher order modes from the 
MLG-like subband, probably because they are masked 
by the much stronger peaks from the BLG-like subband, 
which contains most of the charge density. 
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FIG. 3: Visualizing in-situ band structure changes in TLG with TMF. a-b, Band structure and Fermi surface of 
TLG at zero D (electrostatic potential of each layer equal to zero). The band structure consists of MLG-like and BLG-like 
subbands, with a small band overlap. The bands a and /3 are MLG-like and BLG-like valence bands. The trigonal warping 
effect can be seen in the BLG-like subband. The lattice constant a is 2.46 A. c-d, Band structure and Fermi surface of TLG 
at finite D (for this case, with potential difference between adjacent layers equal to 30 meV). The potential difference induces 
the hybridization between MLG-like and BLG-like subbands and also shifts down in energy the top of the a band, e-h, The 
TMF spectra in TLG at D = 0, 0.18, 0.34, and 0.49 V/nm respectively. As D increases, the a band starts to disappear while 
the /3 band remains visibly unchanged. 



Earlier studies have shown that higher order hopping 
parameters in TLG significantly modify its band struc- 
ture [9l [T8H2T] by introducing subband overlap and trig- 
onal warping in the BLG-like subband (Fig. [3|i). We use 
this full parameter model for the TLG band structure to 
simulate the carrier trajectories and determine the focus- 
ing fields (see supplementary information). The results 
are shown as dashed lines in Fig Although we can 
reproduce the focusing field for the MLG-like subband 
very accurately, we obtain a mismatch in the BLG-like 
subband, similar to those above mentioned in MLG and 
BLG. 

We now focus on the previously unexplored potential 
of TMF as a ballistic electron spectroscopy method to 
investigate in-situ changes in the band structure of a ma- 
terial. One of the remarkable properties of TLG is that 
its band structure can be tuned and controlled by us- 
ing a transverse electric displacement field [22^, D. TMF 
is sensitive to the occupation of each of the TLG sub- 
bands, enabling us to use TMF as a probe of the change 
in the TLG band structure with D. Figures [3^ and b 
show the band structure and Fermi surface of TLG at 
n = —2 X 10^^ cm~^ for the case D = 0. We denote 
the valence bands of the MLG-like and BLG-like sub- 



bands as a and (3 bands, respectively. The application 
of a finite D induces a potential difference between the 
TLG layers, breaking the mirror symmetry and causing 
a hybridization between the MLG-like and BLG-like sub- 
bands. Figures ^ and d show how a finite D results in 
a shift down of the top of the a band . Consequently, 
for a fixed density, the Fermi momentum of the a band 
shrinks with D while that of the (3 band barely changes 
due to its much higher density of states. 

Figures [3^-h show the TMF spectra of TLG at various 
I^'s. We observe a relatively strong focusing peak from 
the a band at = V/nm. However, as D increases, 
the peak starts to shift downward and it eventually dis- 
appears at low density. The disappearance of this peak 
is the result of the top of the a band shifting down in en- 
ergy and leaving the ^ band as a lone contributor to the 
carrier density (Fig. [sj^-d). Therefore, within our density 
range, we end up observing only a single focusing peak, 
from the (3 band, at high D. The onset in the density 
of the focusing peak of the a band allows us to deter- 
mine the potential difference among the TLG layers as a 
function of applied D. As a result, we can estimate the 
effective dielectric constant of TLG which we find to be 
about 3.5 ± 0.2 (see supplementary information). 
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FIG. 4: Temperature dependence of the TMF in MLG and BLG. a-b The TMF spectra as a function of temperature 
from 300 mK to 150 K for MLG and BLG, respectively, at n = —2.8 x 10^^ cm~^. Insets show the amphtudes of the first and 
second modes as a function of temperature (Data taken before current annealing), c transverse magnetic focusing in MLG at 
300 K. The TMF peak as well as the >Jn dependence of the focusing field can be clearly observed (Data taken after current 
annealing) . 



We now look at the temperature dependence of the 
TMF spectra in MLG and BLG. The TMF spectrum is 
affected by temperature, T, at least in two ways: through 
the increase in dephasing (which smoothes the quantum 
interference fluctuations), and through the loss of ballis- 
tic transport due to new scattering channels activated at 
high T. Figures and b show the TMF spectra of MLG 
and BLG, respectively, at n = -2.8 x 10^^ cm~^ from 0.3 
to 150 K. We first concentrate on the fine structure ob- 
served at low T (blue traces) . This structure is the afore- 
mentioned quantum interference between different paths 
on which electrons propagate to the collector |4l [T3HT5] . 
When the temperature-induced broadening of the Fermi 
momentum is on the order of 1/L, electrons become inco- 
herent and the quantum interference is washed out [23], 
resulting in smooth focusing peaks. For our devices, this 
length corresponds to a temperature of about 15 K, which 
is in good agreement with the data. 

In addition, the focusing peaks also decrease as T in- 
creases. The amplitudes of the first and second modes 
are shown in the insets of Fig. and b for MLG and 
BLG, respectively. We observe that the focusing ampli- 
tude in MLG depends linearly on T. A potential scatter- 
ing mechanism includes longitudinal acoustic phonons, 
which give rise to a linear temperature dependence of the 
scattering rate [2 4 . However, we observe a very different 
temperature dependence in BLG. The peak amplitude 



saturates at low temperature and decreases faster than in 
MLG at higher T. A similar temperature dependence of 
the focusing peaks has also been observed in InGaAs/InP 
hetero junctions [25 . Further theoretical work is needed 
to understand the temperature dependence of the focus- 
ing peaks as well as the difference between MLG and 
BLG. 



We end by commenting on the remarkable robustness 
of TMF in graphene. Fig. % shows the TMF fan of MLG 
at room temperature (T = 300 K), where the first mode 
is clearly visible, indicating room temperature ballistic 
transport well into the micron regime. This lower bound 
temperature for the observation of TMF in graphene is at 
least three times higher than the highest temperature at 
which TMF spectrum has been observed in semiconduc- 
tor heterostructures [25^, the main reason probably be- 
ing the lack of remote interfacial phonon scattering [26] 
from hBN. The ability to manipulate ballistic motion 
in graphene at room temperature, coupled with recent 
developments [27] in large area growth of graphene on 
hBN, paves the way towards novel applications based on 
electron-optics. On a more fundamental level, TMF may 
serve as a probe of electron-electron interaction [28H3Q] or 
strain-induced gauge field [3TH33] effects in the electronic 
structure of graphene. 
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Methods 

Figure Ic shows an atomic force microscopy image 
from one of our devices. Our devices are fabricated by 
transferring graphene onto high-quahty hexagonal boron 
nitride [8 . We use oxygen plasma to etch graphene 
flakes into a Hall-bar geometry. Contacts are defined 
by electron-beam lithography and thermal evaporation of 
chromium and gold. The devices are then heat annealed 
in forming gas and subsequently current annealed in vac- 
uum at low temperature [9]. We observe TMF peaks 
both before and after current annealing. The data after 
current annealing have higher quality than before current 
annealing, especially at low density, likely due to reduced 
charge density fluctuations. However, they exhibit simi- 
lar quality at high density. All the data shown here are 
measured after current annealing, except Fig. [2]i-e and 
Fig. [4^-b which was done before current annealing. 

We identify the number of graphene layers by Raman 
spectroscopy and/or quantum Hall measurements. For 
TLG, the quantum Hall measurements reveal that it is 
Bernal-stacked [9j. In addition, we put a top gate onto 
the TLG device, using hBN as a thin dielectric. The com- 
bination of top gate (TG) and bottom gate (BG) allows 
us to control the charge density and the displacement 
field independently. We parameterize the displacement 
field hy D = (CtgVtg + C'BG^BG)/(2eo) where C is the 
capacitive coupling, V is the applied gate voltage rela- 
tive to the charge neutrality point, and eo is the vacuum 
permittivity. 
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